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ABSTRACT 

Here we present the analysis of 3D spectroscopic data of three Blue Compact 
Galaxies (Mrk 324, Mrk 370, and III Zw 102). Each of the more than 22500 
spectra obtained for each galaxy has been fitted by a single gaussian from which 
we have inferred the velocity dispersion (a), the peak intensity (l pea k), and the 
central wavelength (A c ). The analysis shows that the a vs l pea k diagrams look 
remarkably similar to those obtained for giant extragalactic H II regions. They 
all present a supersonic narrow horizontal band that extends across all the range 
of intensities and that result from the massive nuclear star-forming regions of 
every galaxy. The a vs l pea k diagrams present also several inclined bands of 
lower intensity and an even larger a, arising from the large galactic volumes that 
surround the main central emitting knots. Here we also show that the a vs X c 
and A c vs l pea k diagrams, are powerful tools able to unveil the presence of high 
and low mass stellar clusters, and thus allow for the possibility of inferring the 
star formation activity of distant galaxies, even if these are not spatially resolved. 



Subject headings: galaxies: spectroscopy — galaxies: individual (Mrk 370, Mrk 324, 
III Zw 102) -galaxies: dwarfs — galaxies: starbursts 
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Introduction 



High resolution panoramic spectroscopy with good spatial and spectral resolution is 
known to be a powerful tool for studying the kinematics of ionized nebulae as it leads 
to a simultaneous mapping, at seeing limited resolution, of a particular emission line over 
the whole nebula. This however can easily lead to several tens of thousands of spectra, 
making it difficult to issue a detailed and/or a global interpretation of the data. An analysis 
procedure that has proven to be simple and pow erful results from fitt i ng a single gaussian 



to ea ch of the resultant emission line profiles (see IMunoz-Tunon 111994 ; IMunoz-Tunon et al. 



19951 ). A single gaussian regardless of the actual degree of asymmetry or splitting in the 
line profiles. The fit is to conserve the flux of the line profile, and thus lower intensity but 
broader lines would result from the most asymmetric or largely splitted line profiles. From 
the resultant fits one can derive the velocity dispersion (cr), the peak intensity (I pea fc) of the 
fitted lines as well as their central wavelength (A c ). Such a method when applied to giant 
extr agalactic H II region (GHII R) data, leads to two distinct regions in the <j vs ^peak diagram 



see 



Munoz-Tunon et al.lll996l ): A supersonic (cr > chu), relatively narrow horizontal band 



with all possible peak intensity values, and a second region populated by lower intensity 
points presenting even larger supersonic cr values, that crowd along mu l tiple inclined bands. 
For the case of GHIIRs it has been shown (see Munoz-Tunon 1994; IMunoz-Tunon et al. 



19961 : ISabalisck et al.l 1 19951 ; iFuentes-Masip et al.l |2000| ; iTelles et al.l 1200 ll . the last reference 
for the case of H II Galaxies data) that the horizontal band is conformed by data with 
truly gaussian profiles arising mainly from the brightest portions of the nebula, from the 
small regions that enclose groups or clusters of stars. The inclined lower intensity and 
highly supersonic bands on the other hand, emanate from multiple photoionized shells that 
surround the star-forming centers and are likely to result from the stellar mechanical energy 
impact into the ISM, and thus from single gaussian fits to highly asymmetric or strongly 
splitted lines. The analysis of the a vs l pea k diagram for GHIIRs in nearby galaxies (closer 
than 1 Mpc) was then proposed as an excell ent tool to determine their degree of evolution on 
kinematic bases (jMunoz-Tunon et al.lll996l ). Here we extend the analysis to Blue Compact 
Galaxies (BCGs) at even larger distances and show that they experience a similar dynamical 
evolution (ie. their a vs l pea k diagrams look very much like those of GHIIRs). We further 
extend our analysis by looking also at the cr vs A c and A c vs l pea k diagrams, and show them 
as ideal tools to distinctly unveil the presence of high and low mass stellar clusters. 

Here section 2 provides a short description of the observations. Section 3 contains a 
detailed analysis of the data by means of the diagnostic diagrams here proposed and section 
4 gives a summary and a discussion of our findings. 
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The observations 



We use the 6m Big Telescope Alta-azimuthal (BTA) of the Special Astrophysical Ob- 
servatory (SAO) in Nizhnij Arkhyz (Russia), equipped with the multi-mode focal reducer 
SCORPIQjin the scanning Interferometer Fabry- Perot (IFP) observational mode (see lAfanasiev Moiseev 



2005|), for the observations of three BCGs: Mrk 370, III Zw 102, and Mrk 324. SCORPIO 
is provided with a EEV 42-40 CCD array of 2048 x 2048 pixels with an instrumental scale 
of 0.18 " pixel -1 . A set of interference filters centered at the Ha line, covering an interval of 
systemic velocities that range from —200 to 10000 km s -1 , is used at the IFP observational 
mode. The IFP offers two configuration depending on the interference order, in our case at 
the Ha line, the IFP260 (at the interference order 235) and the IFP500 (at the interference 
order 501). Our observation runs were performed using the IFP500 mode, that offers a better 
spectral resolution in the Ha emission line. 

The output from the IFP results on a 3D data cube with x and y being the spatial 
plane and z the wavelength sampling or ethalon step, that in the case of IFP500 mode, it 
translates into 32-40 channels. The original data cube dimensions were (522,522,36) although 
for our analysis we reduced the size, matching the objects size, to (150,150,36). The Free 
Spectral Range (FSR) was 13 A, enough to sample the Ha line, that was scanned with a 
spectral sampling of 0.7 A s=i 32 km s -1 . A binning of 4x4 pixels was applied in the spatial 
directions in order to reduce the time exposure and improve the signal to noise of the final 
data cube, resulting in a final spatial scale of 0.7 " pixel -1 . The instrumental width taken 
from a fitting of the calibration lamp data cube was <Ji ns t = 20 ± 3 km s -1 . Table [1] lists 
a log of the observations together with some data about the galaxies: Column 1 gives the 
name of the object; columns 2 and 3 show the coordinates (from NED0); columns 4 shows 
the linear scale for each galaxy; column 5 lists the total time exposure of the observations; 
column 6 gives the seeing conditions; column 7 shows the systemic velocity of each galaxy 
and column 8 gives the distance in Mpc. Each of our SCORPIO trimmed data cubes contains 
a total of 22500 spectra from an area of ~3 arcmin 2 . The bulk of the data thus requires of 
automatic techniques of analysis to extract the full power of 3D spectroscopy. For each of the 
three galaxies, their Ha emission lines in eve ry pixel were fitted by single gau ssians, to then 



obtain their I pea k, cr, and Ac, as described in iMunoz-Tunon et al.l (119951 . Il996l ). The velocity 



dispersion is given in km s 1 and is corrected for the instrumental and thermal broadening 

2 2 2 "\ 

T obs ~ a inst ~~ a th) 



(a 2 = a 2 bs — <j 2 nst — o- 2 h ). We then used MoisSoft, the software package designed for the 



^dittp: / /www. sao.ru/hq/lsfvo/devices/scorpio/scorpio. html 

2 NASA/IPAC Extragalactic Database is operated by the Jet Propulsion Laboratory, California Institute 
of Technology, under contract with the National Aeronautics and Space Administration(NASA). 
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manipulation of SCORPIO-IFP s pectral line dat a cubes. For a detailed explanation about 
the data reduction procedure, see iMoiseevI ( 20021 ). 



3. Three Dimensional Spectroscopy of BCDs 
3.1. Mrk 324 



Mrk 324 is a galaxy classified as a Blue Compact Dwarf (BCD, see iThuan fc Martin 



198ll ) and is included as a Nucleated BCD in the morphological classification of Cairos 
et al. (2001), as it is basically powered by a luminous central knot. The star formation 
history of the galaxy fits well with an instantaneous star formation law with a Salpeter IMF 
(a = 2.35) and a range of masses between 1< M/M Q <100, in which a central knot of 
4.4 Myr and 1.4 xlO 5 M^ in stars dominates 95% of the emission flux of the galaxy (see 
Martmez-Dekado et aubood ). 



Figure [T] (upper row of panels) shows the results from the single gaussian fit to the full 
Ha SCORPIO data set in the a vs l pea k, cr vs A c , and A c vs l pea k diagram s for Mrk 324. The 
last frame, upper row, shows the Ha image of the nucleated galaxy (see lCairos et al.ll200ll ) 
obtained at the NOT 2.52m telescope of the Roque de Los Muchachos Observatory (ORM) 
in La Palma (Spain). The resemblance of the a vs l pea k diagram with those of GHIIRs is 
remarkable. It presents a horizontal band that limits the values of a to supersonic values (a 
> 20 km s -1 ) and a bunch of low intensity points at even higher a values, that delineate a 
triangular structure above the horizontal band. If one selects from the <j vs ^-peak diagram 
the highest intensity data points, those belonging to the horizontal band (see second row of 
panels in Figured]), these in the a vs A c diagram appear clumped at a given A c (~ 6598.2 
A) what implies that they all arise from a very similar location in the galaxy, that they all 
belong in principle to a single entity. Note however, that this is not exactly coincident with 
the rest wavelength of the galaxy (marked in the A c plots with a red line, see Figured]). 
The selected points in the A c vs l pea k diagram defined also a pointed structure towards high 
intensity values. These points in the emission-line peak intensity map (last frame) display 
their location in the galaxy, which indeed coincides with the major burst of star formation 
in Mrk 324. A major burst with an estimate d mass in stars ~ 4 x 10 5 M ^ and an age of 



about 4.4 Myr (see IMartmez-Delgado M2005I ; iMartmez-Delgado et al.ll2006l ) that causes an H 
II region with a radius of ~ 250 pc. On the other hand, the data belonging to the inclined 
bands with a higher velocity dispersion and a lower intensity in the u vs Ipeafc diagram (see 
third row of panels in Figure [T]) appear well scattered in A c and arise from a large-scale (600 
pc in thickness) expanding rim that surrounds the central emitting knot, as shown in the last 
panel of the third row. The inserts in the last column of Figure [1] display the Ha line and its 
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single gaussian fit at various locations in the galaxy. As in GHIIRs, the data corresponding 
to the horizontal band in the a vs l pea k diagram leads to the highest intensity true gaussian 
profiles, while the profiles fitted to lower intensity and larger o data points, arising from the 
outer expanding rim, are fits to asymmetric and/or double peak line profiles. Figure [T] shows 
the quality of typical line profiles from each of the regions, displaying their intensity, central 
wavelength, and a values. 



3.2. Mrk 370 



This BCD galaxy belongs to the Chained morphological type (see ICairos et al.ll200ll ). 
in which the star formation takes place along a privileged direction. It follows the same 
star formation law as Mrk 324, with 25 young and massive exci ting clusters hoste d over 



an underlying stellar component with ages larger than 5 Gyr, (see lCairos et al.ll2002l ). The 
clusters have an average age of 6.6 Myr with a dispersion of 0.9 Myr, suggesting a coeval 
and global star formation event. For the ones lo cated in the inner nuclear regio n, the stellar 



mass content reaches up to several 10 5 M Q , (see IMartmez-Delgado et al.ll2006l ) 



Figure [2] displays similar diagnostic diagrams for Mrk 3 70. The last frame in the upper 
row of panels shows the Ha image of the chained galaxy (see lCairos et al.ll200ll ) obtained at 
the CAHA 2.2m telescope of the Centro Astronomico Hispano Aleman (CAHA), from Calar 
Alto Observatory (Spain). The ionized gas in Mrk 370 (see the Ha image in Figure[2]) presents 
much more structure than Mrk 324. However, its a vs l pea k diagram looks almost identical 
to that shown in Figure [1] for Mrk 324. The only exception is perhaps the low intensity data 
points with values of a below the supersonic lower limit imposed by the horizontal band (< 
21 km s _1 , see first panel in Figure [2]). The structure of the photoionized gas in Mrk 370 
begins to be disentangled when one selects the data from the horizontal band in the a vs 
^peak diagram (see Figured second row of panels). Such a data set appears clumped at three 
different A c values in the a vs A c diagram. The brightest one sits at the rest wavelength 
of the galaxy (6580. lA). There is a second one at ~6579.4A and a third one between these 
two at ~6579.8A. Two of them are much more apparent in the A c vs I pea k diagram where 
they clearly end up as pointed structures with different intensity (factors of 2) and different 
central wavelengths. The emission-line peak intensity map tracing the selected points (last 
frame in the second row) shows two major bursts of stellar formation and two adjacent much 
less intense, associate d to the massive central s t ellar clusters with masses that range between 
(see 



(8 -50) xlO 4 Mrr 



Martmez-Delgado 1120051 ; IMartmez-Delgado et al.ll2006l ). also resolved 



by ICairos et al.l (120021 ) using the INTEGRAL spectrograph with a single gaussian fit to the 



emission lines. All of these, although near the center of the galaxy, are kinematically resolved 
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as well detached structures. This multiplicity of structures would have been missed if one 
had only use the a vs l pea k diagram as a diagnostic for this galaxy (see first panel in Figure [2]), 
as the brightest points present almost the same supersonic gaussian profiles. These massive 
regions also match the area of the galaxy th at shows the highe st ionization values (see the 



Ha emission-line and the [OlII]/H/3 maps in lCairos et al.ll2002l ). 



The third row of panels in Figure [2] displays the data from the inclined bands with low 
intensity and large a values in the a vs l pea k diagram. This data set, as the corresponding 
one in Mrk 324, presents a large range of values of A c , around the systemic velocity of the 
galaxy. The spread of the data is also noticeable in the A c vs l pea k plane. The bulk of the 
data arises from a large galactic volume shown in the emission-line peak intensity map as an 
expanding rim of gas around the main centers of star formation all over the galaxy. This is 
most intense around the nuclear starburst zone and although less intense, it is still noticeable 
around smaller bursts of stellar formation some 2 kpc away from the nuclear region. 

The final row of panels show the low intensity and low velocity dispersion data from the 
a vs l P eak diagram (see Figure [2]), which appears as vertical well detached low cr structures 
in the a vs A c plane. These structures, given their low values of a, detach well from the bulk 
of the data, most of which presents supersonic values, and due to their largely different A c 
values they also appear as well detached entities along the A scale. Note however that given 
their low intensity values, they are much less apparent in the A c vs l pea k diagram. These 
data points correspond in the emission-line peak intensity map (last panel in the 4 th row 
of Figure [2]) to multiple knots of e mitting gas, or H II regions excited by low mass stellar 



groups (with masses < 10 4 M , see iMartmez-Delgado 1120051 ; IMartmez-Delgado et al.ll2006l ) 



at large distances from the centre of the galaxy. The line profiles arising from these regions 
are also gaussian profiles although with a very low intensity when compared to those arising 
from the central starburst regions (see inserts in the last frame of Figure [2]). 



3.3. Ill Zw 102 



This galaxy is believed to be th e result of an interaction event (see IVorontsov- Velyaminov 



1959: Vorontsov-Vel'Yaminov Il977) and has 



the polar-ring galaxy type (IWhitmore et al. 



jeen considered by some authors to belongs to 



1990). It is morphologycally considered as an 



Extended Blue Compact Galaxy (jCairos et al.ll200ll ). in which 65 star-forming regions were 
identified spreaded over the whole main body of the galaxy, including the arms. The age of 
the knots and their low age d ispersion (6.1 ± 0.6 Myr) suggest a coeval starburst event (see 
Martmez-Dekado et al1l2006h . 
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Figure [3] displays the results for III Zw 102. The last frame in th e upper row of panels, 
shows the Ha image of the extended galaxy (see ICairos et al.l 1200 ll ) obtained at the NOT 
2.52m telescope of the Roque de Los Muchachos Observatory (ORM) in La Palma (Spain). 

The data for III Zw 102 confirms the power of all of the diagnostic diagrams here used. 
The galaxy has a plethora of small H II regions spread along two arms as well as several 
massive centers of stellar formation in the densest nuclear region with m asses between ~ 
10 5 - 10 6 M , fsee lMartmez-Delgadoll2005l : lMartmez-Delgado et al.ll20Q6h . All of these, as 



well as the large galactic volumes undergoing a supersonic expansion leave their signature 
in the diagnostic diagrams. In this way, the horizontal band with supersonic a values at all 
possible intensities noticeable in the a vs l pea k diagram (see Figure 3, 2 nd row of panels) is 
spread over a large range of A c values in the a vs A c diagram and into at least four pointed 
structures (at A c = 6596. 7 A, 6597. 25A, 6598A, and 6600A) of different intensities in the A c 
vs lp e ak diagram. All of these in the emission-line peak intesity map are well resolved as 
giant structures within the galactic nuclear region. On the other hand, the larger a points 
at all peak intensities in the a vs l pea k diagram arise from gas expanding around the nuclear 
zone (see third row of panels in Figure 3). Among this data set is the inclined structure 
between I peafc =1100-1400 which is produced by gas between the four major centers of stellar 
formation in the nuclear region. The broad line structure appears at the center of the galaxy 
with a A c very similar to the systemic velocity of III Zw 102. 

Finally, the lowest a and lowest intensity data points in the a vs l pea k diagram (last row 
of panels in Figure 3) appear, as in the case of Mrk370, as a set of well detached vertical 
stalactites in the a vs A c diagram, hanging from the bulk of the data at very distinct A c values. 
The latter as shown in the last panel are produced by the numerous small H II regions located 
along the galactic arms. As shown by the inserts in the last column of panels, the quality of 
the gaussian fits is totally different in the three regions. The nuclear zone shows supersonic 
perfect gaussians while the surrounding gas presents much broader and asymmetric lines. 
Single gaussians with small o values fit well the data from small centers of star formation. 



4. Summary 



Tridimensional spectroscopy with good spatial and spectral resolution, sampling a par- 
ticular emission line over the whole area covered by an emitting nebula, is known to be the 
most suitable observational technique for analyzing the global kinematics of GHIIRs. In 
particular if a single gaussian fit is carried out over the emission lines in every pixel and from 
those, one infers the velocity dispersion, the central wavelength, and t he peak intensity of 
the lines. The resultant a vs l pea k diagram (IMunoz-Tunon et al.lll996l ) has been shown to 
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be an excellent diagnostic diagram to separate the main broadening mechanisms affecting 
the emission lines, e.g., those that lead to shells and loops generated by the violent action 
from massive stars (see Dyson 1979, Roy et al. 1986) which finally lead to cloud dispersal, 
and the one(s) affecting the regions of massive star formation (either gravity; see Terlevich 
& Melnick 1981, Tenorio-Tagle et al. 1993, or turbulence; see Chu Sz Kennicutt 1994). The 
sequences shown by the various GHIIRs so far analyzed, led the authors also to propose the 
cr-^-peak plot as a tool able to trace the evolutionary status of a GHIIR, from its formation to 
the total dispersal of the ionized gas. 

Here we have confirmed that the data from the inclined supersonic bands in the o vs 
^peak diagrams are caused by the mechanical energy deposited by massive stellar sources. 
These, depending on the age of the stellar clusters and on the density of the surrounding 
gas may affect small or large galactic volumes. However in all cases, the interaction leads to 
very asymmetric line profiles, evident in the a vs l pea k diagram for their highly supersonic a 
values. 

Here we have extended the analysis to Blue Compact Galaxies. We have selected sources 
with a single nuclear starburst (as Mrk 324) as well as galaxies with multiple stellar bursts 
both within their nuclei as well as at large galactocentric distances. In all cases the a vs l pea k 
diagram resembles those from GHIIRs in which the data define a supersonic a (larger than 
Chii ~ 10 km s -1 ) horizontal band with all possible intensity values. Note that the value of 
the limiting a is different in every case. For III Zw 102 the horizontal band peaks at ~ 33 
km s _1 , while for Mrk 324 and Mrk 370 it occurs at 24 and 25 km s l , respectively. 

The a horizontal band, for the case of multiple massive nuclear bursts of stellar forma- 
tion, has been shown to split into different bands in the A c vs l pea k diagram. In this new 
diagram, the largest nuclear regions of massive star formation become evident due to the 
good spectral resolution. And thus even with a poor spatial resolution, or when dealing with 
further away sources, the sampling of the velocity field, with our 3D spectrographs, would 
reveal the exciting sources. Similarly, the a vs A c new diagram picks up the small bursts of 
stellar formation by tracing their lower intensity and slower expansion of their immediate 
surroundings. 

The stalactites detected in the a vs X c diagram in the case of III Zw 102 and Mrk 370 
correspond to the ionized knots located in the outskirts of the galaxies and their spread i n 



a result from the H II regions expansion into a low density medium (see lFranco et al.lll990l ). 



On the whole, the new diagnostic diagrams here presented provide the possibility of 
inferring from 3D spectroscopy, the magnitude of the star formation activity of distant 
galaxies, even if these are not spatially resolved. 
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Fig. 1. — The figure shows in the first row of panels the velocity dispersion a vs peak intensity (the a was corrected for the instrumental and 
thermal broadening), for all the gaussian fits to the individual spectra of Mrk 324. The second and third panels display the corresponding a 
vs central lambda (A c , indicating also the velocity range in km s _1 ) and A c vs peak intensity, respectively. The solid line in these two panels 
indicates the rest velocity of the galaxy. The last panel shows the Ha image of the galaxy and indicates its scale in kpc. The second and 
following rows of panels highlight in colour the various sets of data points selected in the a vs peak intensity diagram and their location in 
the a vs X c and A c vs peak intensity planes. The last panels in every row, with a peak intensity scale, display the locations in the galaxy 
that produce the selected data points. Inserts in these panels correspond to typical gaussian fits (solid red lines) to the data arising from 
different regions (axis correspond to intensity, in counts, vs the wavelength in A). The fitted a values are indicated in every frame. 
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Fig. 2.— Mrk 370. The same as Figured] for Mrk 370. 
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III Zw 102. The same as Figured for III Zw 102. 
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Table 1. Log of the Observations 



Galaxy 


RA(2000) 


DEC(2000) 


Linear Scale(pc/") 


t(s) 


PSF(") 


V H/ (km s- 1 ) 


D(Mpc) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


Mrk370 


02 h 40 m 29 s ?0 


+19°17'50" 


52 


4680 


1.7 


790+ 


10.85 


IIIZwl02 


23 h 20 m 30!l 


+17°13'32" 


110 


6480 


1.8 


1626+ 


22.71 


Mrk324 


23 h 26 m 32!8 


+18°15'59" 


108 


4320 


1.6 


1600+ 


22.43 



+ S ystemic velocity, f r om H I observations, corrected to the Local Group velocity centroid taken 
from lThuan & Martini (|l98lh . 



